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ABSTRACT 

We present deep HST/WFPC2, rest-frame U images of 17 ~ L* quasars at z ~ 1 
and z ~ 2 (V and / bands respectively), designed to explore the host galaxies. We 
fit the images with simple axisymmetric galaxy models, including a point-source, in 
order to separate nuclear and host-galaxy emission. We successfully model all of the 
host galaxies, with luminosities stable to within 0.3 mag. Combining with our earlier 
NICMOS rest-frame optical study of the same sample, we provide the first rest-frame 
U — V colours for a sample of quasar host galaxies. While the optical luminosities 
of their host galaxies indicate that they are drawn purely from the most massive 
(> L*) early-type galaxy population, their colours are systematically bluer than those 
of comparably massive galaxies at the same redshift. The host galaxies of the radio- 
loud quasars (RLQ) in our sample are more luminous than their radio-quiet quasar 
(RQQ) counterparts at each epoch, but have indistinguishable colours, confirming that 
the RLQ's are drawn from only the most massive galaxies (10 11 - 10 12 M Q even at 
z ?s 2), while the RQQ's are slightly less massive (~ 10 11 Mq). This is consistent 
with the well-known anti-correlation between radio-loudness and accretion rate. Using 
simple stellar population "frosting" models we estimate mean star formation rates of 
~ 350 Mq yr^ 1 for the RLQ's and - 100 M Q yr" 1 for the RQQ's at z w 2. By z « 1, 
these rates have fallen to ~ 150 M yr _1 for the RLQ's and ~ 50 M© yr _1 for the 
RQQ's. We conclude that while the host galaxies are extremely massive, they remain 
actively star-forming at, or close to, the epoch of the quasar. 

Key words: galaxies: fundamental parameters - galaxies: active - galaxies: evolution 
- galaxies: high redshift - quasars: general 



1 INTRODUCTION 

Comprehensive studies of low redshift (z < 0.4) quasar 
host galaxies indicate that the most powerful nuclear activ- 
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ity in present-day galaxies is usually as sociated with mas- 

. iDisnev et aT 



sive, bulge-dominated host galaxies (e 



Bahcall et al. 1 19971: [Hooper et al.lll997i; 



McLeod et al 



1995; 



1999: 



McLure et al. 1999; Ridgway et al. 2001; Hamilton et al 
| 2002| : iDunlop et ail [20031; iFlovd et al.l 120041: iGuvon et al 
12009 ; fcanalizo et all l2007l ; TBennert et alJl200Sh . This ties 
in well with the discovery of the black hole - bulge mass re- 
lationship in quiescent galaxies: the host galaxies are exactly 
the type of systems in which one would expect to find the 
very massive black holes that are required to power luminous 
quasars. However, this conclusion is based largely on their 
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Table 1. Summary of the WFPC2 observations: objects at Z f» 1 and 2^2 were observed using the F606W and F814W filters 
respectively (corresponding to approximately rest-frame U in each case). Each quasar is classified as either radio-loud (RLQ) or radio- 
quiet (RQQ) - see Tableland section ^. ll of the text. A single white dwarf star was observed in each filter for PSF calibration. We show 
each target name used in the HST archive, together with its B1950 IAU name. J2000 co-ordinates were obtained from the Digitised Sky 
Survey plates maintained by the Space Telescope Science Institute. Redshifts are derived from NED. Note that the image of SGP5:46 
was contaminated by diffraction spikes from a nearby bright star and had to be excluded from this study. 



Object 


IAU Name 


Type 


2 


J2000 Position 


Integration 


Filter 




(B1950) 






RA (h m s) 


Dec (° ' ") 


time (s) 




SGP5:46 


0049-277 


RQQ 


0.955^ 


00:52:22.8 


-27:30:03 


3600 


F606W 


BVF225 


1301+358 


RQQ 


0.910< 6 ' 


13:04:10.5 


+35:36:51 


3600 


F606W 


BVF247 


1302+361 


RQQ 


0.890( 6 > 


13:05:05.0 


+35:51:21 


3600 


F606W 


BVF262 


1303+360 


RQQ 


0.970( 6 ) 


13:05:30.9 


+35:17:14 


3600 


F606W 


PSF-STAR 


1313+705 


Star 




13:38:59.5 


+70:16:40 


320 


F606W 


PKS0440-00 


0440-003 


RLQ 


0.844^ 


04:42:38.6 


-00:17:43 


3600 


F606W 


PKS0938+f8 


0938+185 


RLQ 


0.940( c > 


09:41:23.2 


+18:21:06 


3600 


F606W 


3C422 


2044-027 


RLQ 


0.942< 6 ) 


20:47:10.4 


-02:36:23 


3600 


F606W 


MC2f f2+172 


2112+172 


RLQ 


0.878W 


21:14:56.7 


+17:29:23 


3600 


F606W 


4C02.54 


2207+020 


RLQ 


0.976< d ) 


22:09:32.8 


+02:18:41 


3600 


F606W 


SGP2:36 


0048-293 


RQQ 


1.773< e > 


00:51:14.3 


-29:05:20 


7800 


F814W 


SGP2:25 


0049-295 


RQQ 


1.869^ 


00:52:07.6 


-29:17:50 


7800 


F814W 


SGP2:fl 


0050-291 


RQQ 


1.976< a ) 


00:52:38.5 


-28:51:13 


7800 


F814W 


SGP3:39 


0053-286 


RQQ 


1.959( e > 


00:55:43.4 


-28:24:10 


7800 


F814W 


SGP4:39 


0056-281 


RQQ 


1.721< e ) 


00:59:08.9 


-27:51:25 


7800 


F814W 


PSF-STAR 


1313+705 


Star 




13:38:59.5 


+70:16:40 


320 


F814W 


PKSf524~f3 


1524-136 


RLQ 


1.687^ 


15:26:59.4 


-13:51:01 


7800 


F814W 


B2 2156+29 


2156+297 


RLQ 


1.753( 6 ) 


21:58:42.0 


+29:59:08 


7800 


F814W 


PKS2204-20 


2204-205 


RLQ 


1.923^) 


22:07:33.9 


-20:38:35 


7800 


F814W 


4C45.51 


2351+456 


RLQ 


1.992(9) 


23:54:22.3 


+45:53:05 


7800 


F814W 



HBovle et all l|l990h : ( ^Hewitt . 



: Burbidgd jl989h: ( c lSchneider et alJ j201Ch: ^ Vcron-Cctty . 
1120011 ): If founlop et alJ d 1989ft ; ^Sticke! fc Kuhj d 1993ft . 



Veronl <1993ft : ^Croom etHl 



luminosities and morphologies in single-band Hubble Space 
Telescope (H ST) imaging (esp ecially at the highest quasar 
luminosities- iFlovd et ajj2004h . The stellar populations and 
masses of quasar host galaxies remain po orly constrained, al - 
though ground-based AT-band imaging llTavlor et all 119961 ) 
and off-nuclear optical spectroscopy (INolan et al.ll200ll ) sug- 
gest the presence (albeit at low level) of an intermediate age 
stellar population in some. It is difficult to disentangle any 
young stellar population in the host galaxy from the glare 
of the active galactic nucleus (AGN), and much of our pic- 
ture of quasar host galaxies is inferred from the properties 
of their inactive, giant red elliptical, counterparts. 

Meanwhile, lower luminosity AGN in the local Universe 
have clear evid e nce f or younger stellar populations (e.g. 
iBoroson fc Okei I1982T ). Seyfert galaxies are frequently 
found with circumnuclear starbursts dHeckman et alJ 



19971 iGu et alJ l200ll: iGonzalez Delgado et al.l l200ll: 



Cid Fernandes et all |2003; iDavies et alj |2007l; iRiffel et al] 



2009). The presence of intermediate age stellar populations 
has also been extensively shown in type 2 (obscured) 
and lo wer luminosity ty p e 1 q u asars llKotilainen fc Ward! 
1994 iRonnback et all Il996l; iBrotherton et all Il999t 



Kauffmann ct al. 2003; Jahnkc et al. l2004allbl ; San chez et aLj 
20041: IVanden Berk et al.1 120061 : Ijahnke et al.l 120071 : 
Jahnke et al.ll2009h . A decrease in mean s tellar age with in- 



creas i ng AGN luminosity is sug gested (e.g. lKauffmann et alJ 
2003; Vanden Berk et aLll2006l ). indicating that we may be 
missing young stellar populations in the most luminous 
nearby quasars, due to the glare of the nucleus. 



At higher red shift we see greater le vels of star forma- 
tion, as expected. iMainieri et al.1 (|201ll ) examined ~ 1800 
obscured (X-ray selected) quasars in the XMM-COSMOS 
survey, with bolometric luminosities > 8x 10 45 erg s _1 . They 
confirm that all dwell in massive (> 10 10 Mq) galaxies, with 
a monotonically increasing AGN fraction with stellar mass. 
They also find strong star formation in the majority, with 
the star forming fraction increasing strongly with redshift 
(62% star forming a t z ~ 1, 71% at z ~ 2, 100% at z ~ 3). 
iTrichas et al.l ([2012) , however, show a decreasing fraction 
of measurable starburst contributions with increasing AGN 
luminosity. 

A key observable is the level of star formation present 
in quasar host galaxies when compared to inactive galax- 
ies of the same mass at the same epoch. In this paper we 
aim to explore the level of star formation in our z ~ 1 and 
z « 2 quasar samples, originally imaged in the rest-frame op- 
tical using HST's Near Infrare d Camera and Mult i-Object 
Spectrometer (NICMOS - see iKukula et all l200ll - here- 
after K01). That paper showed the quasars to be hosted by 
massive (~ 10 11 — 10 12 Mq) elliptical galaxie s, with masses 
consistent with those found at lower redshi ft dMcLure et al.l 
1 19991 : iDunlop et ai1l2003l : IFlovd et al]|2004l - hereafter M99; 
D03; F04 respectively) assuming passively evolving stellar 
populations. Here we present the follow-up rest-frame U- 
band imaging taken with WFPC2, explore the colours and 
morphologies of the host galaxies, and constrain their masses 
and young stellar populations. In particular, we explore their 
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Table 2. Optical and radio properties of the quasars in the 
current study. We show the target names used in the HST 
arc hive, Redshifts, V-ba nd magnitudes, colour excesses E(B — 
V) llSchlegel et al. 1998) are derived from NED. Absolute mag- 
nitudes are confined to the range —23.9 My ^ —25.1 (in rest- 
frame) to ensure that the quasars in each redshift bin are compa- 
rable with each other (as well as with the low-redshift quasars 
in our earlier studies, F04, D03, M99). RQQs have Lsghz < 
10 24 5 W Hz" 1 sr" 1 . RLQs have L 5GH z > 10 24 5 W Hz" 1 sr" 1 
and steep radio spectra. For the estimation of luminosities at 
both radio and optical wavelengths we have assumed a quasar 
spectrum of the form F v oc u~ ' 5 . 

See section r2.1l for full details of sample selection. 



WDJCCTi 


Type 


z 


V 


E(B-V) 


■51.4 








/mag 


1 mag 


/mjy 


Z 


w 1 sample (F606W 'Wide V") 




BVF225 


RQQ 


0.910 


19.6 


0.012 


< 0.15 


BVF247 


RQQ 


0.890 


19.2 


0.014 


< 0.15 


BVF262 


RQQ 


0.970 


19.6 


0.008 


< 0.15 


PKS0440-00 


RLQ 


0.844 


19.2 


0.053 


1443 


PKS0938+18 


RLQ 


0.940 


19.1 


0.030 


432 


3C422 


RLQ 


0.942 


19.5 


0.055 


2117 


MC2112+172 


RLQ 


0.878 


18.7 


0.134 


425 


4C02.54 


RLQ 


0.976 


19.0 


0.045 


790 


z 


« 2 sample (F814W 'Wide /') 




SGP2:36 


RQQ 


1.773 


20.8 


0.014 


< 0.45 


SGP2:25 


RQQ 


1.869 


20.8 


0.016 


< 0.45 


SGP2:11 


RQQ 


1.976 


21.1 


0.017 


< 0.45 


SGP3:39 


RQQ 


1.959 


20.9 


0.028 


< 0.45 


SGP4:39 


RQQ 


1.721 


20.9 


0.022 


< 0.45 


PKS1524-13 


RLQ 


1.687 


21.0 


0.121 


2400 


B2 2156+29 


RLQ 


1.753 


17.5 


0.104 


1223 


PKS2204-20 


RLQ 


1.923 


20.3 


0.049 


700 


4C45.51 


RLQ 


1.992 


20.6 


0.122 


1836 



colours in comparison to the colours of massive galaxies at 
the same epoch. 

The layout of the remainder of the paper is as follows. 
The sample design and observations are discussed in section 
2. We describe the data reduction and analysis in section 3. 
Our results are presented in section 4, and their implications 
discussed in section 5. Our main conclusions are summarised 
in section 6. We assume throughout a cosmology with fl m — 
0.3, A = 0.7 and H = 70 km s _1 Mpc" 1 , and convert 
previous work to this standard where required. 



2 OBSERVATIONS 

Our observing strategy was designed to maximise the de- 
tection of host galaxy light, as in our previous HST studies 
of quasar host galaxies (M99; K01; D03; F04). The precau- 
tions described in those papers (and reiterated below) were 
even more critical to the success of the current observations 
since they target shorter wavelengths than our earlier stud- 
ies. In the rest-frame UV, the quasar spectral energy distri- 
bution (SED) climbs steeply with respect to the host galaxy 
SED, increasing the nuclear-to-host ratio and making accu- 
rate separation of quasar and galaxy light more difficult. The 
observations were performed in HST Cycle 7 with WFPC2, 
as summarised in Table[T] We adopted the Wide Field (WF) 
section of the detector. Despite the higher sampling afforded 



by the PC's smaller (0'.'045) pixels, the O'.'l pixels of the WF 
detectors offered the decisive advantage of better sensitiv- 
ity to low surface brightness emission. The WF3 chip was 
selected for our programme because of its marginally bet- 
ter overall performance relative to the WF2 & WF4 chips 
during the twelve months prior to the observations. 

2.1 Sample design 

Our sample, described in detail by K01, originally consisted 
of 20 quasars divided into two subsamples in the redshift 
ranges 0.83 ^ z ^ 1.0 and 1.67 ^ z ^ 2.01 (referred to 
hereafter as the "z « 1" and "z ~ 2" samples respectively). 
Each subsample contained five radio-loud and five radio- 
quiet objects, and each was matched in terms of its opti- 
cal luminosity-redshift distribution. The sample was defined 
to fall in the luminosity range —24 ^ Mv ^ —25, allow- 
ing them to be compared directly to each other and to the 
quasars in our lower- redshift samples (M99; D03; F04). How- 
ever, the results presented by K01 were at odds with this 
criterion, with one z ~ 2 quasar shining close to Mv = —27, 
and the z ~ 2 RQQ's being significantly lower in total opti- 
cal luminosity than their RLQ counterparts. This appears to 
be partly due to errors in the (highly heterogeneous) photo- 
metric measurements used to define the sample (which pre- 
dated SDSS and 2QZ), as well as uncertain Tf-corrections 
from observed to rest-frame V at signifi cant redshift s. 

The RQQs were derived from the iBovle et al.l 
UVX survey ("SGP" objects) or iMarshall et all 
survey ("BVF" objects), and confirmed as radio 



1990) 



1984) 



Hz _1 sr _1 ) by the 



quiet 

Very Large Ar- 



ray (VLA) FIRSTS and NVSS3 surveys. RLQs were se- 
lected from the Veron-Cetty & Veron (1993) quasar cat- 
alogue, with reported 5-GHz radio luminosities L^ghz > 
10 25,5 W Hz -1 sr -1 and steep radio spectra to ensure that 
their intrinsic radio luminosities are not boosted by relativis- 
tic beaming. The NICMOS observations of two quasars in 
the original K01 sample (the RQQ "SGP2:47" at z = 0.830 
and the RLQ "1148+56W1" at z = 1.782) were affected by 
technical problems, and so these were not re-observed with 
WFPC2. We also note that the WFPC2 image of SGP5:46 
was contaminated by diffraction spikes from a nearby bright 
star and had to be excluded from the present study. 



2.1.1 Present Sample 

The current sample therefore consists of eight objects at 
z » 1 and nine at z w 2, for which the radio and optical 
properties are listed in Table [2] The updated redshift versus 
rest-frame optical luminosity distribution for our sample is 
shown on the left of Fig. [T] using the HST rest-frame V pho- 
tometry presented in K01. The QSO's span t he knee of the 
quasa r luminosity function at each redshift l|Croom et al.l 
l2004h . and the radio quiet subsamples are still reasonably 
well matched to the optical luminosities of our z « 0.2 and 
0.4 samples (M99, D03, F04). The shrinking of, and biases 
inherent within, our sample affect our original aim of explor- 
ing the differences between the RLQ and RQQ populations. 



1 http://sundog.stsci.edu/ 

2 http://www.cv.nrao.edu/nvss/ 
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Figure 1. Left: Absolute V magnitude versus redshift for the quasars in our host galaxy studies out to z ?s 2. Filled circles represent 
radio-quiet quasars, while open circles represent radio-loud quasars. Objects in the present WFPC2 study at z 1 & 2 arc shown using 
large symbols, compared to our lower redshift objects (F04, D03, M99). Our sample spans the knee of the quasar lumin osity function at 
each redshift: the dotted line indicates M v for the 2QZ quasar luminosity function at each redshift llCroom et al.|[2~004l ) . See section \1. II 
for notes on the sample selection. Right: Illustration showing generic spectra for a quasar nucleus (upper panel) and an early-type 
galaxy (lower panel), with our approximate rest-frame bandpasses marked. We have tailored our filter selection to target the rest-frame 
U and V-band, thus sampling the SED of the host on either side of the break feature at 4000A. While the Mg II quasar emission line is 
admitted, we avoid prominent galaxy emission lines (section l2.2H . 



We focus instead on providing a benchmark colour for all 
quasar host galaxies at these redshifts. We note that the 
conclusion drawn by KOI that RQQ hosts are less luminous 
than RLQ hosts, is potentially biased at z ~ 2 by the a pos- 
teori luminosity distribution of the quasars themselves (see 
Fig- HI- This is discussed further in sections 14.41 and 15.51 



2.2 Choice of filters 

We chose filters that maximise sensitivity to the underlying 
starlight in the wavelength range of interest, while avoiding 
strong galaxy emission lines that might positively bias our 
measurement of the host galaxy luminosity. Our NICMOS 
observations sampled the rest-frame V-band continuum at 
each redshift, avoiding the strong emission lines [Oiii]A5007 
and Ha using the F110M (w J) filter at z rts 1 and F165M 
(» H) at z « 2. With WFPC2 we targetted the objects' 
rest-frame near-ultraviolet (UV) continuum, in the region 
~ 2300 — 3400A, corresponding closely to U band. We use 
WFPC2 F300W as our reference rest-frame U band through- 
out this work. 

For the z « 1 sample we have used the F606W (Wide- 
V) filter, and for z w 2 the F814W (Wide I) filter in order 
to maximise throughput. We have avoided strong galaxy 
emission lines such as [O ll]A3727. The absence of strong 
galaxy emission lines in the observed bandpass is impor- 
tant as they can render the images sensitive to any large 
scale regions of emission-line gas which might be present. 
These are directly associated with the nuclear activity and 
can mask or confuse our picture of the underlying starlight. 
However, the quasar broad emission line Mg nA 2796 is in- 
cluded by using such wide filters (see Fig. [T] right panel). 
The equiv alent width of Mg H in the SDSS composite quasar 
spectrum jVanden Berk et al.ll200l1 ) is ~ 32A, and will con- 
tribute a slight (up to 0.03 mag) positive bias to our nu- 
clear luminosities. However, since Mg II arises very close to 
the nucleus and such high ionization lines are only found 



in absorption in galaxies, our host galaxy luminosities are 
unaffected. 



2.3 Extinction and iC-corrections 

For the estimation of luminosities at both radio and op- 
tical wavelengths throughout this work, we have assumed 
a quasar spectrum of the form F v ~ v~ ' 5 . it-corrections 
are applied to convert to rest-frame magnitudes (see sec- 
tion 2]3J), where we revisit this assumption. We have con- 
verted the original counts from the quasar images of KOI to 
AB mags for comparison in this paper. For galactic extinc- 
ti on corrections we use colour excess values from the map 
of lSchlegel et al.l (|l99ct ) and a diffuse interstellar dust model 
(Rv ~ 3.1). All magnitudes are AB (so that a flat energy 
distribution in F v corresponds to zero colour). 



2.4 Saturation and integration times 

Integration times were calculated to allow the WFPC2 im- 
ages to reach a similar depth, and to detect the host galaxies 
out to a similar radius, as in the previous NICMOS observa- 
tions. To do this we made the null hypothesis that the host 
galaxies are passively evolving. Any star formation would 
thus serve to boost the galaxies' luminosities, making them 
easier to detect. We therefore observed each of the z ~ 1 
quasars for two orbits through the F606W filter, and each 
of the z ~ 2 quasars for four orbits through the F814W 
filter. Each HST orbit was sufficient to include three long 
exposures (600 or 700 s), plus several snapshots. The lat- 
ter provided a good characterisation of the bright nuclear 
point source in the event that this region of the image was 
saturated in the longer integrations. This resulted in a total 
on-source integration time of 3600 s for the z » 1 objects, 
and 7800 s for those atzsi2 (see Table EE). 
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2.5 Point Spread Function 

Experience from our previous WFPC2 & NICMOS imag- 
ing studies of quasar host galaxies (M99; KOI; D03, F04) 
has demonstrated that it is essential to obtain a deep, high- 
dynamic range stellar Point Spread Function (PSF) through 
each filter in order to accurately characterise the critical 
outer regions of the quasar. Synthetic TINYTI1V0 PSFs, 
whilst an excellent match to the inner ~ 1 arcsec, fail to 
adequately reproduce the complex and variable outer struc- 
ture of the PSF wings. Archive PSFs are rarely precisely cen- 
tred on the same region of the detector as the target quasars 
(thus changing the structure of the outer region of the PSF), 
and do not always have the requisite dynamic range. In any 
case the form of the WFPC2 PSF changes systematically 
with time and so measurements made as close to the epoch 
of the observations as possible are highly desirable. For this 
study we therefore adopted the same strategy as used in our 
previous HST programs. An entire orbit was devoted to ob- 
taining high-dynamic range stellar PSFs through each of our 
two chosen filters. The star selected (GRW+70D5824) is a 
white dwarf, with a flat spectral shape over the wavelengths 
of interest which is very similar to that of a typical quasar, 
thus ensuring the best possible match to the nuclear PSF 
of each target source. Into the centre of this deep, empirical 
stellar PSF we splice a TINYTIM realisation of the theoret- 
ical PSF at the location of the centre of the quasar on the 
detector. The quality of fit is strongly dependent on the fo- 
cus term in the TINYTIM model, which affects the amount 
of light in the central peak. Focus changes significantly dur- 
ing an orbit due to thermal "breathing" of the telescope 
structural. We tested each observation with a range of fo- 
cus values by performing a \ 2 minimisation on a fit of the 
PSF to the central region of each quasar (best fits accept- 
able at 3cr level). In each case the adopted PSF had a focus 
close to the value determined for the telescope at the time 
of observation. 

3 DATA REDUCTION AND ANALYSIS 

3.1 Data reduction 

Data reduction was carried out using the HST/WFPC2 
pipeline and standard procedures in iraf. An additional 
step was required for the images of the PSF star and those 
quasars for which the nucleus had saturated in the longer ex- 
posures (specifically, all nine quasars in the sample at z ~ 1 
plus SGP2:25 from the z ~ 2 sample). The saturated regions 
of these images were replaced by the corresponding regions 
(suitably scaled) from the unsaturated snapshot image of 
the same object. 

3.2 Two-dimensional modelling 

In order to separate the galaxy and nuclear components of 
the quasar image we used the 2-D modelling software devel- 
oped for our previous host-galaxy studies (se e M99; F04 for 
details) and tested and described in detail bv lMcLure et alj 

3 http:/ /www. stsci.edu/software/tinytim/tinytim. html 

4 http: / / www.stsci.edu /hst / observatory/focus / 



|2000l) ; iFlovdl (|2005h : iFlovd et al.l l|200Sl ). The code uses 
downhill simplex \ 2 minimization to match a synthetic 
galaxy plus central point source to the HST image, with 
nuclear luminosity and galaxy effective radius (i? e ), surface 
brightness (/i e ), axial ratio (a/b) and position angle (9) as 
free parameters. The model galaxy's surface brightness pro- 
file can be set to either an exponentially decaying disc, or an 
r 1//4 de Vaucouleurs law characteristic of an elliptical galaxy. 
The model is convolved with an appropriate PSF, described 
in section [2751 above. 

Accurate characterisation of the uncertainties in the 
data is critical when fitting across such high dynamic range. 
The transition from PSF sampling error (see F04 for details 
of sampling error treatment) to Poissonian error is deter- 
mined by finding the radius at which the sampling error 
drops to 1.5 times the mean background Poisson noise. The 
noise on the central pixel was fixed by the Poisson statistic 
for that pixel. 

We fit the host galaxy in each image with both an ellip- 
tical galaxy (r 1//4 -law profile) and an exponential disc galaxy 
model. The difference in the \ 2 values of the two fits is used 
to classify the host as disc or bulge-dominated. Sersic, and 
bulge-disc decomposition models are not fitted (c.f. F04) as 
we are unable to constrain the additional parameters. Each 
fit is performed out to ~ 3" for the z ~ 1 objects and ~ 2" 
at z ~ 2 (determined from the radial profiles). As discussed 
further in section [4] we find that an Elliptical host galaxy 
morphology provides the best fit in all cases. 

Note that we adopt the circularly-averaged effective ra- 
dius, R e (that is the mean radius of the half-light elliptical 
isophote) throughout this work. We convert KOl's half-light 
semi-major axis lengths to this standard. This ensures that 
the size measurement is not positively biased by high eccen- 
tricities. 



3.3 Fixed R e models 

In addition to the conventional modelling, we also fitted each 
object using a series of fixed-radius models (at R e — 2, 5 
and 10 kpc). This technique was adopted by KOI for many 
of their z ~ 2 objects, since they were unable to recover 
reliable measures of the R e at this redshift. We chose to 
apply the technique to all the objects in the sample in order 
to determine the stability of the host-galaxy luminosities 
recovered by the modelling, and in order to compare host- 
galaxy models of identical radii at different wavelengths. 



4 RESULTS 

The best fit galaxy and nuclear parameters derived from the 
2D modelling procedure are listed in tables [3] (for the z w 1 
sample) and [5] (for the z ~ 2 sample). We have successfully 
recovered stable (±0.3 mag) host-galaxy fluxes for all seven- 
teen quasars studied. The stability of the host and nuclear 
fluxes are shown by the \ 2 contours in Fig. [2] Host-galaxy 
V — J (z w 1) and I — H(z~2) colour data derived from the 
WFPC2 (this paper) and NICMOS (KOI) images are pre- 
sented in tables[S]and[S]respectively. The WFPC2 images of 
the quasars, their best-fit models, modelling residuals and 
radial profiles are presented along with brief notes on each 
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Figure 2. Ax 2 contours (at 1, 2 and 3<r for the host galaxy - solid; and nucleus - dashed) in the m,jj st ~ m Nuc plane. All magnitudes 
are AB. The dotted lines indicate constant total flux. The host-galaxy fluxes are constrained to within ±0.3 mag, and the nuclear fluxes 
generally somewhat better. Those that deviate (MC2U2+172, B2 2156+29, PKS2204-20 and 4C45.51) show significant circumnuclear 
residual flux (Figs. lA"71lA14|[A16IIA17ll . that is confused with the nucleus, thus weakening our constraint on its total flux (see section [5.111 . 
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Table 3. Best fit models for the quasars at 2 rj 1 (all Elliptical). Columns are as follows: object name; reduced-x 2 value for the best 
fit model; A\ 2 between the elliptical and disk-morphology model; number of degrees of freedom in fit, u; effective radius, R e , of best 
fitting galaxy model in kpc; surface brightness of the host at the effective radius, fj, e , in units of Vab mag.arcsec -2 ; integrated apparent 
magnitudes of the nucleus and the host galaxy in VAB-band; the ratio of integrated nuclear and host galaxy luminosities; position angle 
of the host (in degrees East of North); the axial ratio of the host. All fluxes are presented as apparent AB magnitudes, uncorrected for 
galactic extinction. 



Object 


Y 2 


Ax 2 


V 


R € 
/kpc 


He 

/mag.arcsec - 


T/nuc 
AB 

2 /mag 


T/host 

/mag 


Ln/Lh 


PA 

1° 


a/6 












Radio-Quiet Q 


uasars 










BVF225 


1.122 


134 


1526 


5.5 ±0.2 


23.9 ±0.3 


18.5 ±0.1 


21.2 ± 0.1 


12.6 


176 


1.24 


BVF247 


1.143 


97 


1395 


4.5 ±0.2 


24.5 ±0.5 


20.4 ±0.1 


22.3 ±0.1 


6.1 


155 


1.25 


BVF262 


1.077 


82 


1512 


3.0 ±0.1 


23.2 ±0.5 


20.0 ±0.1 


22.0 ±0.1 


5.9 


4 


1.33 












Radio-Loud Q 


uasars 










PKS0440 


1.097 


55 


1797 


11.1 ± 0.5 25.5 ±0.4 


18.9 ±0.1 


21.3 ±0.1 


9.2 


47 


1.12 


PKS0938 


1.106 


45 


2195 


3.8 ±0.2 


23.3 ±0.6 


19.3 ±0.1 


21.5 ± 0.1 


7.4 


154 


1.20 


3C422 


1.073 


24 


2812 


2.7 ±0.1 


22.2 ±0.4 


19.2 ±0.1 


21.2 ± 0.1 


6.2 


87 


1.57 


MC2112 


1.040 


3 


1806 


1.8 ±0.1 


20.5 ±0.1 


19.9 ±0.1 


20.2 ± 0.1 


1.4 


102 


1.74 


4C02.54 


1.070 


31 


2360 


6.3 ±0.1 


24.2 ±0.1 


18.8 ±0.1 


21.4 ±0.1 


10.2 


111 


1.42 



Table 4. Best fit models for the quasars at z 2 (all Elliptical). Columns are as follows: object name; reduced-x 2 value for the best 
fit model; A\ 2 between the elliptical and disk-morphology model; number of degrees of freedom in fit, v\ effective radius, R e , of best 
fitting galaxy model in kpc; surface brightness of the host at the effective radius, £t e , in units of 7ab mag.arcsec -2 ; integrated apparent 
magnitudes of the nucleus and the host galaxy in /^B-band; the ratio of integrated nuclear and host-galaxy luminosities; position angle 
of the host (in degrees East of North); the axial ratio of the host. All fluxes are presented as apparent AB magnitudes, uncorrected for 
galactic extinction. 
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v 2 , 
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Re 
/kpc 


He 

/mag.arcsec -2 


rime 
AB 

/mag 


rhost 
AB 

/mag 


Ln / Lh 


PA 
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a/6 










Radio-Quiet Quasars 










SGP2:36 


1.244 


115 


1164 


1.6 ±0.1 


22.3 ± 0.4 


21.5 ±0.1 


22.5 ±0.1 


2.5 


136 


1.29 


SGP2:25 


1.295 


61 


1106 


1.8 ±0.1 


23.2 ± 0.1 


20.8 ±0.1 


23.1 ±0.2 


8.3 


98 


1.21 


SGP2:11 


1.169 


5 


1101 


4.0 ±0.8 


25.2 ± 2.0 


20.6 ±0.1 


23.4 ±0.3 


12.3 


168 


1.15 


SGP3:39 


1.325 


51 


1189 


2.2 ±0.2 


22.8 ±0.3 


20.6 ±0.1 


22.5 ±0.1 


5.7 


107 


1.13 


SGP4:39 


1.031 


2 


1041 


1.5 ± 1.9 


22.1 ± 2.1 


20.7 ±0.1 


22.5 ±0.2 


5.1 


81 


1.09 












Radio-Loud Quasars 










PKS1524 


1.198 


39 


1121 


2.2 ±0.2 


22.0 ±0.6 


19.6 ±0.1 


21.5 ±0.1 


5.5 


163 


1.31 


B2 2156 


1.335 


45 


1110 


3.6 ±0.1 


22.3 ±0.1 


19.4 ±0.1 


20.7 ±0.1 


3.4 


39 


1.10 


PKS2204 


1.300 


25 


981 


1.6 ±0.1 


21.3 ±0.2 


20.8 ±0.1 


21.6 ±0.1 


2.0 


137 


1.48 


4C45.51 


1.767 


50 


1119 


2.2 ±0.1 


22.0 ±0.2 


20.6 ±0.1 


21.5 ±0.1 


2.2 


78 


1.13 



object in Appendix A (Figs. lA"TrlA17[) . Finally, the fixed R e 
model fits are presented in Appendix B (tables |BT1 fc IB2| ) . 

The fits to the z ~ 1 objects are formally acceptable 
at the 4ct level. In general, the z ~ 2 objects have best fit 
models that are not formally acceptable, and the resulting 
parameter estimates are therefore uncertain. However, we 
find that the host-galaxy luminosity is very stable to changes 
in the fitted parameters. This can be seen both from the x 2 
contours (Fig. [2} , and from the results of the fixed R e fits 
(tables [BT1fc [B2)) . 

4.1 Host-galaxy morphologies 

As mentioned above, in all cases a bulge-dominated host 
galaxy provides a marginally better fit (up to 3a) to the 
data than a disc-dominated host (see tables [3] and |4}. We 
note that all of the z ~ 1 objects were robustly determined 
to be elliptical in the optical by K01. For this reason, as well 
as our marginal morphological preference, we show the best- 
fit elliptical galaxy model in each case. We note, however, 



that the smooth r 1//4 -law profile does not form as good a 
fit to the data as in our previous studies (M99; K01; D03; 
F04). All the images exhibit significant "bumps" in their 
radial profiles within 2" of the nucleus (see Figs. IA1TtA17[) . 
These bumps are directly traceable to numerous features in 
the images which may be nearby companions, mergers, or 
knots of star formation within the host galaxies themselves. 
We have masked these features out from the fit in order 
to model just the underlying host-galaxy light distribution. 
Aside from these caveats, in all cases the best-fit disc and 
bulge models do give consistent (within 3cr) values for the 
total flux from the host galaxy, giving confidence that this 
parameter has been determined robustly. 

4.2 Fixed R e fits 

The results of the fixed R e modelling are presented in Ta- 
bles [Bl] and [B2] along with the best- fit models obtained by 
K01 for reference. Note that K01 presented semi-major axis 
scale lengths, whereas we use effective radii, R e (mean ra- 
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Table 5. z es 1 host-galaxy and nuclear V — J colour, in AB 
magnitudes, after correction for galactic extinction. 



Object 


V-J 


V-J) 




(Nuc) 


(Host) 


Radio-Quiet Quasars 


BVF225 


-0.1 ±0.3 


0.5 ±0.3 


BVF247 


-0.4 ±0.3 


2.8 ±0.3 


BVF262 


0.2 ±0.3 


1.5 ±0.3 


Radio-Loud Quasars 


PKS0440 


-0.3 ±0.3 


1.8 ±0.3 


PKS0938 


-1.2 ±0.3 


1.3 ±0.3 


3C422 


0.6 ±0.3 


2.2 ±0.3 


MC2112 


0.0 ±0.3 


1.1 ±0.3 


4C02.54 


0.5 ±0.3 


1.3 ±0.3 



dius of the half-light elliptical isophote - see section 13.21) . 
We note that several of KOl's host-galaxy scale lengths ap- 
peared artificially large due to their high eccentricity. We 
have converted scale lengths to R e in the tables. 

Generally the UV and optical morphologies of the z ~ 1 
host galaxies are indistinguishable (best fit morphology, R e , 
axial ratio and position angle). However BVF247, 3C 422 
and MC2112 exhibit compact UV host galaxies compared 
to the optical (and noticeable circumnuclear residuals). At 
z ~ 2 the UV host galaxy detections are quite compact com- 
pared to the optical. In each case, adjusting R e over a rea- 
sonable range makes little or no difference to the luminosity 
of the best-fit host galaxy. Furthermore, adopting the fixed- 
R e model that best matches each NICMOS (K01) model in 
place of our own preferred model results in no significant 
changes to the results. All of the following conclusions are 
thus robust to the uncertainty in the true size of the host 
galaxy. 

4.3 Observed fluxes and colours 

The observed optical-IR fluxes and colours for the quasar nu- 
clei and host galaxies, after correction for galactic extinction, 
are shown in tables [S] and O NIR fluxes are from the models 
of K01 (converted to AB), while optical fluxes come from 
our modelling results above. The mean extinction-corrected 
apparent host-galaxy colours for the quasar sample are (with 
associated standard errors): 

(V- J) zxl = 1.9 ±0.3 

{I-H) z „ 2 = 1.1 ±0.3 

For the nuclei they are: 

(V- J) z(a i = -0.1 ±0.2 

{I-H) zx2 = 0.4 ±0.2 

4.4 A"-corrections and intrinsic luminosities 

The observational bandpasses were chosen to coincide 
closely with rest-frame U and V for our WFPC2 and NIC- 
MOS observations, respectively. We calculate the U and V 
band luminosities by /("-correcting the observed fluxes into 



Table 6. z ~ 2 host-galaxy and nuclear I — H colour, in AB 
magnitudes, after correction for galactic extinction. 



Object 


I — H 


I — H 




(Nuc) 


(Host) 



Radio-Quiet Quasars 



SGP2:36 


0.2 ±0.3 


1.5 ±0.3 


SGP2:25 


-0.1 ±0.3 


1.9 ±0.4 


SGP2:11 


0.4 ±0.3 


1.4 ±0.4 


SGP3:39 


-0.3 ±0.3 


1.4 ±0.3 


SGP4:39 


0.5 ±0.3 


-0.5 ±0.4 


Radio-Loud Quasars 


PKS1524 


0.1 ±0.3 


0.6 ±0.3 


B2 2156 


0.0 ±0.3 


1.4 ±0.3 


PKS2204 


0.9 ±0.3 


-0.4 ±0.3 


4C45.51 


1.7 ±0.3 


2.2 ±0.3 



these emitted bands (see lHogg et al.ll2002l ) using the band- 
pass curves for each filter. For both the host galaxies and the 
QSO's themselves we assume power-law SED's, F v cx v~ a , 
or equivalently F\ oc \^ , where /3 = a — 2. For the host galax- 
ies we adopt a po wer law, Phost = —0.39, a ppropriate for a 
distant red galaxy (jvan Dokkum et al . 2006). For the QSO's 
we have assumed /3qso = —1.5 (a — 0.5) as defined earlier 
(section |2,3[) . and consistent with our average measured nu- 
clear colours, although it is clear that there is actually a large 
range of variation from object to objec t, with a recent study 
provi ding a range a = 0.55 ± 0.42 l|Kennefick fc Bursickl 
120081 ). Since our observed bandpasses coincide closely with 
the reference emitted bands, the exact slope of the SED for 
each quasar does not strongly affect our luminosity calcula- 
tion, it just affects the distribution of flux across the filter 
bandpass. In fact, varying the index over a very wide range 
of realistic values (— 1 < /3 host < 1; — 2 < /3qso < — 1) yields 
only ±0.05 mag. uncertainty on the host-galaxy luminosity 
and ±0.03 mag. uncertainty in the quasar nuclear luminos- 
ity. 

The absolute U and V band magnitudes for the quasar 
host galaxies and nuclei are presented in Table[7] As in K01, 
we find that the hosts of radio-loud quasars are more lumi- 
nous than those of the radio-quiet quasars, and strongly so 
at z ~ 2 (see section [5T)l . However, at z ~ 2 this is likely to 
be due to the optical luminosity biases in the RLQ (optically 
brighter) and RQQ (dimmer) subsamples. As discussed in 
section |2~T1 the sample was originally designed to match the 
optical luminosity and redshift distributions of the RL and 
RQ subsamples, but was based on old and heterogeneous 
photometric measurements, with large A"-corrections from 
observed to rest-frame V. The results of K01 (their Fig. 6 
and our Fig. [TJ show that the subsamples are not matched 
in luminosity. This is further discussed in section IS~5l 

The host galaxies' average rest-frame U — V colours 
(with associated standard error of the mean) are as follows: 

(U-V) x mi = 1.6 ±0.3 

(U- V) z „2 = 0.9 ±0.3 

In spite of the optical and radio luminosity differences be- 
tween the RL and RQ samples, the colours of their host 
galaxies are indistinguishable in the present data. 
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Table 7. Absolute magnitudes and rest-frame U — V colours for the quasar host galaxies and nuclei at z ~ 1 and 2. U magnitudes are 
F300W and V are Subaru V (all AB, corrected for galactic extinction). 



Ubject 


lVLxj (IN UC j 


My (Host) 


Mv(Nuc) 


My (Host) 


TT T/ 

U — V 


TT T/ 

U — V 












(IN UC ) 


/TT j \ 

(rlOSt ) 






Radio-Quiet Quasars 






BVF225 


-24.7 ±0.1 


-21.9 ±0.1 


-24.6 ±0.3 


-22.4 ±0.3 


-0.1 ±0.3 


0.48 ±0.3 


Dvr Z4t ( 


OO *7 _1_ f| 1 
— Zz. f ± U.l 


-20.7 ±0.1 


-22.3 ±0.3 


-23.5 ± 0.3 


— U.4 it U.O 


o ei _|_ n q 


r> V r zbz 


OQ Q _l_ n 1 

— 23.3 ±0.1 


-21.3 ±0.1 


-23.4 ±0.3 


-22.8 ± 0.3 


0.2 ± 0.3 


1.49 ± 0.3 






Radio-Loud Quasars 






.rlvbU44U 


o A o _J_ n 1 
— 24.2 ± U. 1 


-21.7 ±0.1 


-23.9 ±0.3 


-23.4 ±0.3 


—0.3 ± 0.3 


1.75 ± 0.3 




o /l n _l_ n 1 
— 24. U ± U.l 


-21.8 ±0.1 


-22.8 ±0.3 


-23.1 ± 0.3 


— 1.2 ± 0.3 


L.oZ ± U.O 


3C422 


O A O _l_ A 1 

— 24.2 ± U. 1 


-22.2 ±0.1 


-24.8 ±0.3 


-24.3 ± 0.3 


0.6 ± 0.3 


2.15 ± 0.3 




oq a _l_ n 1 
— 23. ± U.l 


-23.2 ±0.1 


-23.6 ±0.3 


-24.3 ± 0.3 


0.0 ± 0.3 


1.14 ± 0.3 


4L/UZ.D4 


— Z4.0 ± U. 1 


-22.1 ±0.1 


-25.1 ±0.3 


-23.4 ±0.3 


U.O ± U.O 


1 QO I n Q 






Radio-Quiet Quasars 






SGP2:36 


-23.1 ±0.1 


-22.1 ±0.1 


-23.2 ±0.3 


-23.4 ±0.3 


0.1 ±0.3 


1.25 ±0.3 


SGP2:25 


-23.9 ±0.1 


-21.7 ±0.2 


-23.7 ±0.3 


-23.4 ±0.3 


-0.1 ±0.3 


1.72 ±0.4 


SGP2:11 


-24.2 ±0.1 


-21.6 ±0.3 


-24.5 ±0.3 


-22.8 ± 0.3 


0.3 ±0.3 


1.20 ±0.4 


SGP3:39 


-24.2 ±0.1 


-22.5 ±0.1 


-23.9 ±0.3 


-23.6 ± 0.3 


-0.3 ±0.3 


1.19 ±0.3 


SGP4:39 


-23.8 ±0.1 


-22.1 ±0.2 


-24.3 ±0.3 


-21.4 ±0.3 


0.4 ±0.3 


-0.68 ±0.4 






Radio-Loud Quasars 






PKS1524 


-25.0 ±0.1 


-23.2 ±0.1 


-25.0 ±0.3 


-23.6 ± 0.3 


0.0 ±0.3 


0.44 ±0.3 


B2 2156 


-25.3 ±0.1 


-24.0 ±0.1 


-25.3 ±0.3 


-25.2 ± 0.3 


-0.1 ±0.3 


1.21 ±0.3 


PKS2204 


-24.0 ±0.1 


-23.3 ±0.1 


-24.8 ±0.3 


-22.7 ±0.3 


0.8 ±0.3 


-0.62 ±0.3 


4C45.51 


-24.4 ±0.1 


-23.7 ±0.1 


-26.0 ±0.3 


-25.7 ±0.3 


1.6 ±0.3 


1.98 ±0.3 



5 DISCUSSION 

5.1 Are we detecting host-galaxy light? 

Our observations were designed to detect the UV light from 
even a mature stellar population, and it is clear directly 
from the images themselves (see Figs. IA1I - IA17[) that in 
general we have an excess of extended UV flux (compared 
to the PSF). However, in three cases at z ~ 1, and every 
case at z ~ 2, the rest-frame (7-band host-galaxy effective 
radius is significantly smaller than the rest-frame V-band 
equivalent reported by K01 (see tables [rBll fc IB2 (l . Are we 
truly detecting the host-galaxy light, or are we being biased 
by the strength of the AGN, or perhaps by reflected AGN 
emission? 

We find that the host-galaxy and nuclear luminosities 
are quite stable to changes in the effective radius of the host 
galaxy. The fixed R e models (tables IB 1 1 & IB2|I show a sig- 
nificant increase in x 2 as we move away from the preferred 
radius to larger or smaller values, while the host-galaxy lu- 
minosity remains stable. However, the stability of the host 
and nuclear luminosities is demonstrated most clearly by the 
X 2 contours in m,Host — rn^ UCl (Fig. [2}. These show that the 
host-galaxy flux is constrained to within ±0.1 to ±0.3 mag. 
The nuclear luminosity is generally somewhat better con- 
strained (±0.1 to ±0.2 mag). We also note that the host 
galaxies generally have similar axial ratios and position an- 
gles in each band (tables [BT1fc[B2)l. 

In almost all cases, the nuclear U — V colours are close 
to those expected for a typical QSO SED. Only 4C45.51, 
PKS2204 and PKS0938 (all radio-loud) strongly deviate 
from the expected flat nuclear colour. We ran PSF-only fits 
to all the objects, and found a significantly higher \ 2 for the 
best PSF-only fits in every case. Finally we note that while 



our host galaxies are compact, they are still well resolved, 
and have comparable sizes to the ma ssive (> 10" Mq ) 
galaxies at 1 < z < 3 in CANDELS ijBruce et al.ll2012u . 
This is discussed further in section 15.41 

We are therefore confident that we are constraining the 
total host-galaxy luminosity in all cases. However, our R e 
and /x e constraints are sometimes weak, and the 7? 1 ^ 4 -law 
model is only marginally preferred over an exponential disk 
at z as 2. There is clearly UV-bright substructure in the 
host galaxies that diminishes the quality of a smooth model 
profile fit to the data. 



5.2 Residual features 

We note the significant residual features in MC2112±172 
(Fig. [A7|l B2 2156±129 (Fig. fATo) . and PKS2204-20 
(Fig. IA16|) . The former two are objects with compact U 
band hosts compared to their V band hosts (PKS2204— 20 
had a poorly constrained host-galaxy size in K01). These 
residuals are not consistent with PSF features, being more 
extended than the Airy ring and off-centre with respect to 
the nucleus. Our pixel scale corresponds to ~ 800 pc at the 
distance of our sources, and the feature seen in MC2112±172 
is up to w 10 kpc across. This may represent an asymmet- 
ric disk component in the host galaxy or a foreground or 
merging companion. B2 2156+129 shows a compact clump 
of excess flux slightly off-centre (to the northwest), up to 
« 4 kpc across. PKS2204— 20 shows a slight "bridge" of ma- 
terial across the nucleus, linking up with the companion to 
the northwest. 4C45.51 shows a very compact slightly off- 
centre (to the north) feature that is ~ 2 kpc across, although 
this may be PSF residual. These features explain the unusual 
contours for these objects' fits (Fig. [2}: the relatively poor 
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Figure 3. Host-galaxy observed colour with redshift for the RQQs (filled circles) and RLQs (open circles) compared to COSMOS galaxies 
(density plot) and SSP models (lines). Magnitudes are all extinction-corrected AB in the COSMOS filters. Left: V — J host colours for 
the qu asars at z fa 1. Right: I — H host colours for the quasars at z fa 2. SSP colour evolution curves are shown for a lBruzual fc Charlod 
(2003) passively evolving population (red line), plus "frosting" models with from 0.01% (orange) up to a maximum (purple) of 5.12% 
of t he galaxy mass in ongoing star formation (in octaves). All SSP models assume a redshift of formation zj = 5, solar metallicity 
and lChabrien i2003T l IMF - see section [53] All COSMOS galaxies with acceptable best- fit photometric redshifts out to z = 5 are shown. 
Quasar host galaxies are bluer than passive and have similar colours to the full range of COSMOS galaxies at similar redshifts. 




Figure 4. Host-galaxy observed colour with optical host luminosity (K01) for the RQQs (filled circles) and RLQs (open circles) compared 
to COSMOS galaxies within 3<r of the mean quasar redshift (density plot - see section [5.31 1. Also show are the best linear fit (dotted 
line) and median colour at each luminosity (dashed line). Magnitudes are all AB in the COSMOS filters. Left: V — J host colours for 
the quasars at z fa 1. Right: I — H host colours for the quasars at z fa 2. Quasar host galaxies are bluer than expected for galaxies of 
the same redshift and luminosity - all fall below the median COSMOS colour (two borderline at 2 fa 2). 



constraint on the nuclear flux is due to confusion between 
the nucleus and the contaminant. Although these are the 
most spectacular examples, all of our images show some ex- 
cess flux above the smooth model, and the majority (all at 
z fa 2 and three at z fa 1) show a clear non-PSF, centrally 
concentrated component. 

5.3 Host-galaxy luminosity and colour 

We now compare our quasar host-galaxy colours and lumi- 
nosities with the colours and luminosities of typical galaxies 
at the same redshifts (Figs. [3] and [4]). We adopt the data of 
the COSMOS galaxy survey, since it is sufficiently deep and 
wide to provide a statistical sample of galaxies at the same 
redshift s as our quasars. We use the C OSMOS galaxy cat- 
alogue l|Capak et al.ll2007l ; lilbert et al.ll2009l ) for redshifts, 
V, J and I band data, and the COS MOS _g-band cata - 
logue (P. Capak, private communication; [Bielbv et al.ll2012j ) 
to compare observed colours with our host-galaxy sample. 



We adopted the best fit photometric redshift, but removed 
all objects for which the 1-a relative uncertainty in the red- 
shift was greater than 10%, or for which the quality of the 
redshift fit was poor (x 2 /v > 2). We have converted our lu- 
minosities and fluxes to the equivalent COSMOS bandpasses 
for ease of comparison. 

Fig. [3] shows the observed galaxy colour against red- 
shift for all such COSMOS galaxies out to z — 5 (density 
plot), and for our quasar host galaxies (circles). We also 
show the colour evolution of several synthetic simple stellar 
population (SSP) m odels computed using the s tellar spectral 
synthesis models of lBruzual fc Charlotl |2003l ). We show an 
instantaneous burst of star formation at a redshift Zf = 5, 
wi th solar metalli city and the initial mass function (IMF) 
of lChabried (|2003l ) that includes low mass / sub-stellar ob- 
jects yielding more accurate mass-to-light ratio determina- 
tions. The red line in each case shows the evolution of a pas- 
sively evolving stellar population. "Frosting" models, based 
on the same underlying old stellar population are shown 
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Figure 5. Mean rest-frame U — V colour, and V band luminosity, versus redshift, for the quasar host galaxies of RQQ's (filled circles) and 
RLQ's (open circles) at Z w 1 and z « 2. Sample means (by type and by redshift bin) are shown with standard errors. All magnitudes are 
AB, V -band is the COSMOS V-band and U-band is F300W. Luminosity and colour evolution curves are shown for a lBruzual fc Charlotl 
1120031 1 passively evolving SSP (red line), plus "frosting" models with 0.01 (orange), 0.02, 0.04, 0.08, 0.16, 0.32, 0.64, 1.28, 2.56 and 5.12% 
(purple) of the galaxy mass in ongoing star formation. Luminosities are shown for s tellar mas ses of 10 11 (solid lines) and 10 12 M (dotted 
lines). All SSP models assume a redshift of formation zj = 5, solar metallicity and lChabrierl j2003T ) IMF. Quasar host-galaxy colours are 
consistent with ~ 0.1% star formation at 2 fa 1 and ~ 1% at z fa 2. The colours of the RLQ and RQQ host galaxies are indistinguishable 
at each epoch, in spite of the large mass differences present in the sample. 



in colour, with 0.01, 0.02, 0.04, 0.08, 0.16, 0.32, 0.64, 1.28, 
2.56 and 5.12% of the galaxy mass undergoing star forma- 
tion. The young stellar population is modelled with a single 
burst of 0.01 Gyr age, and our estimates of the active star 
forming fractions are therefore lower limits (see section I5~4")) . 

In Fig. [4l we plot the observed galaxy colour against 
luminosity. Here, comparison samples of COSMOS galaxies 
are selected to be within 3<r of the mean redshift of our 
samples; that is 0.780 < z < 1.058 at * « 1 and 1.505 < 
2 < 2.196 at 2 « 2. This cut is in addition to the quality 
cut described above. The COSMOS data illustrate the well 
known correlation between galaxy luminosity and colour. A 
blue cloud and red sequence are clearly apparent at z w 1. A 
hint of the red sequence is still visible at z « 2, although at 
this redshift the most massive galaxies are no longer purely 
red. We also plot the best linear fit to the COSMOS colour- 
luminosity data (dotted lines) and the mean colour in 11 
equal sized luminosity bins (dashed line). 

The quasar host galaxies are extremely bright, sampling 
the upper end of the galaxy luminosity function (~ 1 — 4L*). 
The majority of such luminous COSMOS galaxies are found 
to be red at z ~ 1. At z ~ 2 there is a less clear-cut colour- 
magnitude division due to low numbers. However, all of our 
quasar host galaxies are bluer than the median galaxy colour 
for COSMOS galaxies of the same luminosity and redshift 
(Fig. [4]— dashed line). The probability of all 17 quasar host 
galaxies being bluer than the median, if they are drawn at 
random from the galaxy distribution at similar redshift and 
luminosity, is 0.0008%. Allowing for the two objects that 
have error bars that admit redder-than-median colours, the 
probability (15/17) is 0.12%. This result is robust to changes 
in the host-galaxy scale length, and to matching R e to the 
best fit models in K01. We therefore conclude that the host 
galaxies of quasars are not drawn at random from the mas- 
sive galaxy population, but are systematically bluer than 
the average galaxy of the same mass and redshift. 



5.4 Host-galaxy luminosities, masses and 
star-forming fractions 

Fig. [5] shows the host-galaxy rest- frame U — V colours and 
V band luminosities. For comparison we show the evolving 
luminosities and colours of the same set of SSP models as 
discussed in the previous section. Luminosity evolution is 
shown for stellar masses of 10 11 (solid lines) and 10 12 M Q 
(dotted lines). From their colours, the quasar host galaxies 
at 2 ~ 2 are 0.5 — 5% star-forming by mass. By z » 1 this 
has fallen to 0.05 - 0.5%. 

The resulting host-galaxy masses are somewhat lower 
than would be estimated from the simple passive curves (e.g. 
K01), since there is significant contribution to the optical 
luminosity from the young stellar population. This is seen 
most clearly in Fig. [5] The blue colour of the host galaxies 
implies a shallower slope to the luminosity evolution, and 
thus a significantly lower mass than would be inferred as- 
suming a purely passive model. Note that the mass difference 
is more extreme than one might naively imagine for the few 
percent "frosting" fractions, since the young stellar popula- 
tion contributes strongly to the optical luminosity relative 
to its low mass. Host-galaxy stellar mass ranges at z w 2 
are: 

• 1 x 10 11 - 1.3 x 10 12 M for the RLQ's. 

• 3 x 10 10 - 3 x 10 11 M for the RQQ's. 

By 2 w 1 this has increased to: 

• 4 x 10 11 - 1.3 x 10 12 M for the RLQ's. 

• 2 x 10 11 - 6 x 10 11 M for the RQQ's. 

There is roughly an order of magnitude scatter in the 
estimated mass of each subsample thanks to the range of 
colours seen. This spread in mass dominates over the ex- 
pected ~ 0.3 — 0.6 dex sy stematic error in the mass estimates 
from our SSP modelling (|Gallazzi fc Bellll2009l ; IConrov et all 
2009). The masses are shown against effective radius in 
Fig. El togeth er with the massiv e 1 < z < 3 CANDELS 
galaxies from iBruce et al. (2012). As we remarked in sec- 
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tion l5.1l above. our host-galaxy models are quite consistent 
with this sample of massive galaxies at similar redshifts. 
Our z ~ 1 host galaxy effective radii are, on average, a 
fa ctor of 2 larger than those for the lower luminosity AGN 
of ISanchez et all i|2004h in the same redshift range, but with 
significant overlap in the populations. 

We have modelled the young component with a O.OlGyr 
old stellar population, and our colours and luminosities sug- 
gest star formation rates ~ 350 Mq yr _1 for the RLQ's 
and ~ 100 Mq yr" 1 for the RQQ's at z « 2. By z « 1, 
these rates have fallen to ~ 150 M© yr " 1 for the RLQ's and 
~ 50 M yr _1 for the RQQ's. The iKennicutt j"r] (ll998T l 
relation for the UV luminosity of star formation gives simi- 
lar estimates. Modelling the young component with a range 
of ages (0.001-0.1 Gyr) gave consistent results for the star- 
formation rate. 

The mean mass of the RLQ's grows ~ 10% between 
z » 1 and z ~ 2. However, the upper limit on the masses 
of the RLQ host galaxies remains reassuringly unchanged 
at each redshift, as this is already 10 12 M© at z ~ 2. This 
gives the illusion that these objects are passively evolving, 
although it is clear from their colours at each epoch that 
this is not the case. It is simply that at each epoch, the 
RLQ's exist in the most massive galaxies, and these have 
already assembled as much stellar mass as the most massive 
galaxies known by z ~ 2. There is, however, an increase 
in the host-galaxy masses of the RQQ's observed between 
z ~ 2 and z ~ 1. An apparent mass difference between the 
RLQ and RQQ host galaxies remains at each redshift, while 
the star-forming fractions appear indistinguishable (but see 
section [531 below). At both epochs, there is plenty of time 
for the host galaxies to evolve onto the red sequence by 
the present day, but to avoid the production of extremely 
massive galaxies in the present day Universe, it is clear that 
the star-formation seen in these objects is not sustainable. 



5.5 Radio-Loudness and Accretion Rate 



As mentioned earlier (sections 12 . 1.11 l4~4")l . the RLQ and RQQ 
samples were initially selected to have identical observed op- 
tical luminosities. However, the results of K01 showed the 
actual rest-frame optical luminosities to be lower for the 
RQQ's at z « 2, and the host galaxy luminosities to be 
significantly lower at both redshifts. This is consistent with 
the well-known anti-correlation between radio-lou dness and 
accretion rate, seen in lower luminosity AGN (seelXu et al.l 
1 19991 ; lHoll2002l ; iPanessa et aHl2007l : ISikora et al.ll2007h . Our 
results demonstrate that this anti-correlation holds even for 
the most optically luminous RLQs and RQQs. RLQs (shin- 
ing at lower Eddington ratios) will have higher black hole 
masses for a given optical-UV luminosity than their RQQ 
cousins (which have higher Eddington ratios). Through the 
well-known scaling relations between black-hole and galaxy 
bulge mass, we will therefore also observe a higher mass 
population of host galaxies for the RLQs. 

Given this difference in accretion rate between radio- 
loud and radio-quiet quasars, the similarity in the colours 
of their host galaxies (and hence in their specific star for- 
mation rates) is perhaps all the more remarkable. While the 
present data are insufficient to draw strong conclusions, we 
note that this is consistent with optical quasar activity be- 
ing associated a particular phase in the evolution of massive 




Mass/M 

Figure 6. Mean galaxy size versus mean stellar mass for our host- 
galaxy samples (large circles for z as 1; triangles for z ^ 2; open 
for RLQ; closed for RQQ), co mpared to th e > 10 11 Mq, 1 < 
z < 3 CANDELS galaxies of iBruce et alj [I2012T ). Uncertainties 
shown are standard errors on the mean. While compact, our host- 
galaxy detections are quite consistent with the typical sizes of 
massive galaxies at the same redshift. 



galaxies, for example shortly after a starburst, or shortly 
after truncation of ongoing star formation. 



5.6 Interpretation and future observations 

Quasar host galaxies at z = 1 — 2 appear to be exceptional, 
both in mass (as at lower redshift) and in their colour for 
their mass (which has been a contentious matter for the past 
decade at lower redshifts). Our results therefore fit in well 
with a growing consensus that quasar and AGN activity is 
strongly associated with ongoing star formation, although it 
is too early to establish cause and effect. 

Studies of the host galaxies of l ower luminosity AGN 
and quasars at low redshift [z w 0.2 - Ijahnke et alj 2004alV 
intermediate redshift (0.5 < z < 1.1 - ISanchez et al.ll2004h 
and high redshift (1.8 < z < 2.75 - Ijahnke et alj|2004bh 
present a consistent picture, which this study confirms and 
extends to higher AGN luminosity. All of the above stud- 
ies find a mix of morphological types for the host galaxies 
(as expected at lower AGN luminosities). However, where 
the host galaxies are found to be elliptical, their colours are 
consistently found to be bluer-than-average for inactive el- 
liptical galaxies. 

Our results are consistent with those of lMainieri et all 

l|201ll ) who showed that the star forming fraction in power- 
ful X-ray selected quasars (Lboi > 8 x 10 45 erg s _1 ) increases 
strongly with redshift, with 62% star forming at z ~ 1, 71% 
at z ~ 2, 100% at z ~ 3. However, the quasars in this study 
would appear to be more strongly star forming than those 
at lower redshift, where such activity has been detected. By 
stacki ng radio images of COSMOS X-ray selected AGN at 
2 < l. lPierce et al.l (|201lT ) showed that the sources detected 
at 24/im also have radio emission, and are consistent with 
star formation rates of up to 40 M© yr _1 . Evidence of star 
formation in quasars is also given by the recent detection 
of large amounts (~ 10 10 Mq) of cold molecular gas in the 
majority of IR-ultraluminous QSO's, and the strong corre- 
lation of CO luminosity with FIR lumin osity for all QSO's, 
IRQSO's and ULIRG's jXia et al.ll20ll ). 
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We support the findings of IShi et al.1 l|2009h who pro- 
pose that type 1 quasars reside in a distinct population of 
galaxies that exhibits elliptical morphology but harbours a 
significant fraction of intermediate-age stars a nd is expe- 
riencing intense circumnuclear star formation. ICenl (|2012T ) 
propose a model for coevolution of galaxies and supermas- 
sive black holes whereby a central starburst precedes the 
quasar phase, which is fuelled by remnants of the starburst. 
A major prediction is that luminous quasar host galaxies 
are a special population of early-type galaxies located in the 
green valley of the galaxy colour-magnitude diagram. Our 
data are broadly consistent with this model, although our 
data suggest a rather bluer colour for the host galaxies than 
would be expected if they lie exclusively in the green val- 
ley We certainly cannot reject it due to possible sources of 
UV luminosity bias (UV-bright clumps in the host galaxy, 
and a possible contribution from scattered AGN light). The 
star-formation levels detected here, as a fraction of galaxy 
m ass, are an order of magnitude greater than those detected 
bv lNolan et ail l|200lT ) in z ~ 0.2 quasar host galaxies (af- 
ter correcting for their use of a 0.1 Gyr old young stellar 
population). It will be interesting to revisit the low redshift 
samples in the UV to pick out the locations of the young 
stellar populations and compare them to the ti d al fea tures 
observed in the deep imaging of iBennert et al.1 (|2008l ). Fi- 
nally, our esti mated star-formati on rates are consistent with 
those found bv lPage et all | |2012D , although our z ~ 2 RLQ's 
would appear to be slightly more strongly star forming than 
their uppe r limit s (~ 350 Mq as oppose d to < 300 Mq. 
iPage et al.l < |2012f ) and lTrichas et al.1 (|2012h posit the lack of 
the very highest star-formation rates in quasar host galax- 
ies as evidence for AGN feedback. Our similar range of star 
formation rates, and the similarity in colour of the RLQ 
and RQQ host galaxies support this picture. It remains pos- 
sible that the star formation triggers the quasar activity, 
but in that case we might expect to see even stronger star- 
formation in the host galaxies. 



6 CONCLUSIONS 

The main result of this study is that UV-luminous host 
galaxies are present in our entire sample of radio-loud and 
radio-quiet quasars at z ~ 1 and z ~ 2. Our fitted host- 
galaxy luminosities are stable to within ±0.3 mag. and our 
conclusions are robust to changes in the preferred UV host- 
galaxy size. They occupy elliptical galaxies in the upper re- 
gions of the optical luminosity function (1 — 4L*), but have 
bluer than median colours for galaxies of the same redshift 
and luminosity. They are strongly star forming at z ~ 2 
(~ 1% by mass), with weaker (~ 0.1%) star- formation in 
somewhat more massive galaxies by z ~ 1. The RLQ's 
and RQQ's are indistinguishable in terms of colour, but the 
RLQ's are more massive at each epoch, indicative of their 
lower accretion rates (see section 15. 5[) . The most massive 
RLQ's have stellar masses of ~ 10 12 Mg, even at z ~ 2. 

The probability of obtaining galaxies as UV-bright as 
these quasar hosts from a sample of quiescent galaxies of 
the same optical luminosity and redshift is low (< 0.12%), 
and we reject the hypothesis that quasar host galaxies are 
drawn at random from the massive galaxy population. We 
conclude that while the host galaxies are extremely massive, 



they remain actively star-forming (or very recently so) at the 
epoch of the quasar. The relationship between AGN activity 
and star-formation remains unknown, but it seems clear that 
optically luminous AGN activity is associated with recent or 
ongoing star-formation in massive galaxies, whereas radio 
AGN activity is associated purely with the most massive 
galaxies. 

The host galaxies are more compact than low-redshift 
analogues, but are consistent with massive quiescent galax- 
ies at same epoch. There is clear evidence for dumpiness in 
the UV flux. We reiterate that the models are designed to 
fit the smooth bulk of the starlight. However, we have noted 
that the radial profiles of the host galaxies presented here 
are not as smooth as those in the optical. The UV "bumps" 
correspond to compact regions of increased UV luminosity - 
perhaps close companion objects or knots of star formation 
in the host galaxy. These were not apparent in the NICMOS 
rest-frame V imaging. The UV dumpiness, and compact UV 
morphology observed in some objects suggest that higher 
resolution imaging and IFU spectroscopy may reveal com- 
pact regions of star formation in the host galaxy, and hint 
at the possibility of discovering circumnuclear starbursts in 
some particularly UV-compact objects. Our results there- 
fore support the growing consensus that quasar and AGN 
activity is strongly associated with recent star formation 
(possibly its truncation - see section [53]) , and suggest that 
this activity tends to be located in the central regions of the 
galaxy. 
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APPENDIX A: IMAGES, MODELS AND 
NOTES ON EACH OBJECT 

In this section we present the following for each object (left- 
to-right and moving downwards): A 10" square "postage 
stamp" image centred on the quasar; The best fit convolved 
model; The model-subtracted residuals; The radial profile. 
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Figure Al. The RQQ BVF225 at z = 0.910. A prominent nu- 
cleus and an interacting companion are the most noticeable fea- 
tures of this object. The field contains a number of other small 
galaxies greater than 5 arcsec from the quasar. It is strongly nu- 
clear dominated in the rest-frame optical ([imic/^host] v = 7.6) 
and in the rest-frame U ([L nuc / L^ ost ]i/ = 12.6). Significant PSF- 
related circumnuclear flux is visible in the residuals. 

In the radial profile, the best-fit model is shown as a solid 
line, with the PSF component indicated by a dotted line. 
Data points are azimuthal averages with la error bars. Note 
the significant deviations from smooth profiles in the major- 
ity of cases. These bumps are traceable to numerous features 
in the image which may be nearby companions or mergers, 
or knots of star formation within the host galaxies them- 
selves. 

This paper has been typeset from a TjnjX/ I^TjjX file prepared 
by the author. 
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Figure A2. The RQQ BVF247 at z = 0.890. A number of can- 
didate companions are visible, most noticeably a small, possibly 
interacting object some 2" to the North. Whilst the host is quite 
prominent in rest-frame V , it is heavily nuclear dominated in the 
UV: (L nuc /L host ) v = 0.3; (L nuc / L host ) v = 6.1. There is some 
excess circumnuclear flux. 
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Figure A3. The RQQ BVF262 at z = 0.970. Nuclear dom- 
ination is more extreme in U as expected: [L nuc / L^ oat ] y = 
1.8; [L nuc /L hoBt ]u = 5.9. 
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Figure A4. The RLQ PKS0440-00 at z = 0.844. A high- 
powered radio quasar, with a Gigahertz-peaked spectrum. We 
note a number of very compact candidate companions in the 
field, but there are no obvious mergers or interactions. The ob- 
ject is more nuclear-dominated at U as expected: [inuc/ihost] V = 
1.4; [inuc/ihostlc; = 9-2 



Figure A6. The RLQ 3C422 at z = 0.942. A typical compact 
steep-spectrum object. There is one obvious, comparably large 
companion 6" to the south, as well as a number of smaller objects 
dotted around the field. The quasar is nuclear- dominated at U : 
[L nuc /L host ] v = 1.4; [L nuc /L host ]u = 6.2. Some PSF-related 
circumnuclear flux is visible in the residuals 
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Figure A5. The RLQ PKS0938+18 at z = 0.940. We see the 
expected increase in nuclear domination by a factor of ~ 3 from 
V to U: [L nuc /L hoBt ] v = 2.3; [L nuc / L hoBt ] v = 7.4. There is a 
loose collection of small candidate companions out to a radius of 
around 20". In addition, there is some very low surface brightness 
"fuzz" , 3" NW of the nucleus, which is left over in the modelling 
residuals and may indicate a merger. 
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Figure A7. The RLQ MC2112+172 at z = 0.878. An elon- 
gated object in our WFPC2 image, suggesting a significant on- 
going merger event, possibly with a small object visible to the 
NNE. There are a number of small candidate companions in the 
field. The residuals show significant (apparently non-PSF) ex- 
cess flux west of the nucleus. The host galaxy is prominent, with 
a slightly stronger nuclear component in U: [L nuc / L^ ost ] y = 
0.5; [Lnuc/ihost](7 = 1-4. 
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Figure A8. The RLQ 4C02.54 at z = 0.976. In our image there 
is clear evidence for a merging companion 3" to the northwest, 
in addition to a number of other galaxies in close proximity to 
the quasar. If the large object 5" to the southeast is found to 
be at the same redshift as the quasar, then it too is a large el- 
liptical galaxy. We find a higher level of nuclear dominance in 
rest-frame U: [L nU c/£host] V" = 4.8; [L nU c/£host][/ = 10.2. There 
are significant residuals to the immediate north of the nucleus. 
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Figure A9. The RQQ SGP2:36 at z = 1.773. K01 found a rel- 
atively weak nucleus in V: [Imuc/Lhost] V = 0.8. In U we again 
find a relatively weak nucleus, with [L nuc / Lh os t](7 = 2.5. A large 
number of small objects are visible in the field. 




-4 -2 2 4 -4 -2 2 4 




-4 -2 2 4 



Figure A10. The RQQ SGP2:25 at z = 1.869. This quasar is 
quite close on the sky to the quasar SGP2:27 at z = 1.93. K01 
found a moderate nuclear component; [L nuc / L^ ost ]y = 1.3. We 
find a much stronger nuclear component; [L nuc / L^ ost }u = 8.3. A 
large number of faint objects are observed in the field. 
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Figure All. The RQQ SGP2:11 at z = 1.976. We were un- 
able to constrain the size or surface brightness of the host, but 
the total host-galaxy luminosity is well constrained. The quasar 
is nuclear-dominated at both wavelengths: [L nuc / L^ OBt ]y = 
4.7; inuc/^hoat]c/ ~ 12.3. There are some circumnuclear resid- 
uals, but these are compact. 
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Figure A12. The RQQ SGP3:39 at z = 1.959. SGP3:39 lies close 
on the sky to the lower redshift quasar SGP3:35 (z = 1.498). 
The quasar has a fairly weak nuclear component in V, but much 
stronger in U: [L nuc /L boat ] v = 1.2; [L nuc / L host ]u = 5.7. 
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Figure A14. The RLQ PKS1524-13 at z = 1.687. We see a large 
number of foreground / companion objects in the field, particu- 
larly to the south. The nucleus dominates at both wavelengths: 
[inuc/^host] V = 3.2; [L nuc / Lhostlu = 5-5- 
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Figure A13. The RQQ SGP4:39 at z = 1.721. K01 found it to be 
a heavily nuclear-dominated quasar in V ([L nuc / L^ ost ]y = 12.4). 
In U the host-galaxy size and surface brightness are poorly con- 
strained due to the strong nucleus ([L nuc / L^ oat ]u ~ 5). However, 
the total host-galaxy magnitude is well constrained. The residu- 
als are remarkably clean. The apparent tail linking the quasar to 
the companion 3" west is probably diffraction spike residual that 
happens to coincide with the companion. 
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Figure A15. T he RLQ B2 2156+2 9 at z = 1.753. A number of 
nearby galaxies {Thomas et alll 995) are visible in the image. The 
nucleus dominates in U: [L nuc / L hoBt ] v = 0.9; [L nuc / L hoBt ] v = 
3.4. There is a significant circumnuclear residual, likely a PSF 
artefact. 
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Figure A16. The RLQ PKS2204-20 at Z = 1.923. There are two 
quasars with this PKS identifier, this one being the more distant 
one. The other, at z = 1.62, is out of our field of view. Unusually, 
the nuclear domination in rest-frame U appears to be less than 
that in the optical: [L nuc / L host ] v = 3.2; [L nuc /L host ]u = 2.0. 
There is a notable companion 2" NW of the nucleus that ap- 
pears to be interacting with the host galaxy. 
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Figure A17. The RLQ 4C45.51 at z = 1.992. [L nuc /L h ost] v = 
1.4; [i n uc/ihost](7 = 2.2. We see diffraction spikes from a nearby 
star at the top (southeast) of the image. There are four notice- 
able foreground / companion objects within 10" and the residuals 
show a compact (2 pixel) artefact just north of the nucleus that 
may be attributable to the PSF. 
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Table Bl. Fixed R e results for the z ~ 1 sample. The best-fit 
model from Table \3\ is presented, together with models at fixed 
i? e = 2, 5 and 10 kpc, and the best fit (J-band) model from 
KOI. Note that the KOI sizes have been updated for the current 
cosmology and converted to R e . 



Table B2. Fixed R e results for the z ss 2 sample. The best-fit 
model from Table [4] is presented, together with models at fixed 
Re = 2, 5 and 10 kpc, and the best fit (ff-band) model from 
KOI. Note that the KOI sizes have been updated for the current 
cosmology and converted to R e . 



Object 




Rc 




t/Nuc 
AB 


T/Host 
AB 


PA 


a/b 


BVF225 


1.122 


5.5 


23.9 


18.5 


21.2 


176 


1.24 




1.189 


2.0 


21.2 


18.5 


20.9 


4 


1.20 




1.134 


5.0 


23.5 


18.4 


21.2 


178 


1.23 




1.130 


10.0 


25.0 


18.4 


21.3 


171 


1.29 


J-band 




5.0 


24.4 


18.5 


20.7 


33 


3.20 


BVF247 


1.143 


4.5 


24.5 


20.4 


22.3 


155 


1.25 




1.194 


2.0 


22.1 


20.7 


21.9 


154 


1.26 




1.189 


5.0 


24.3 


20.6 


22.0 


156 


1.25 




1.209 


10.0 


25.7 


20.6 


22.0 


157 


1.26 


J-band 




9.2 


23.4 


20.7 


19.5 


3 


1.20 


BVF262 


1.077 


3.0 


23.2 


20.0 


22.0 


1 


1.33 




1.164 


2.0 


22.3 


19.9 


22.1 


3 


1.44 




1.153 


5.0 


24.5 


19.9 


22.3 


7 


1.52 




1.155 


10.0 


26.0 


19.9 


22.3 


10 


1.61 


J-band 




3.1 


22.2 


19.8 


20.4 


5. 


1.36 


PKS0440 


1.097 


11.1 


25.5 


18.9 


21.3 


17 


1.12 




1.160 


2.0 


21.2 


19.2 


20.9 


43 


1.05 




1.133 


5.0 


23.5 


19.1 


21.2 


48 


1.08 




1.137 


10.0 


24.9 


19.1 


21.1 


47 


1.10 


J-band 




7.3 


23.3 


19.1 


19.4 


2 


1.64 


PKS0938 


1.106 


3.8 


23.3 


19.3 


21.5 


154 


1.20 




1.110 


2.0 


21.5 


19.3 


21.3 


135 


1.22 




1.108 


5.0 


23.8 


19.3 


21.6 


1 10 


1.25 




1.118 


10.0 


25.3 


19.3 


21.6 


1 12 


1.28 


J-band 




3.0 


21.7 


20.4 


20.1 


136 


1.33 


3C422 


1.073 


2.7 


22.2 


19.2 


21.2 


87 


1.57 




1.076 


2.0 


21.2 


19.2 


21.0 


88 


1.58 




1.079 


5.0 


23.8 


19.1 


21.5 


89 


1.73 




1.090 


10.0 


25.6 


19.0 


21.8 


91 


2.05 


J-band 




11.9 


23.5 


18.5 


18.9 


139 


1.33 


MC2112 


1.040 


1.8 


20.5 


19.9 


20.2 


102 


1.74 




1.040 


2.0 


20.4 


20.2 


20.2 


101 


1.65 




1.148 


5.0 


22.9 


20.0 


20.7 


102 


1.69 




1.265 


10.0 


24.6 


20.0 


20.8 


102 


1.81 


J-band 




15.8 


23.7 


19.6 


18.8 


168. 


1.02 


4C02.54 


1.070 


6.3 


24.2 


18.8 


21.4 


Ill 


1.42 




1.102 


2.0 


21.1 


19.0 


20.9 


112 


1.26 




1.071 


5.0 


23.5 


18.9 


21.3 


112 


1.36 




1.072 


10.0 


25.1 


18.8 


21.4 


111 


1.50 


J-band 




2.4 


22.2 


18.2 


19.9 


16. 


4.08 



APPENDIX B: 
RESULTS 



FIXED EFFECTIVE RADIUS 



Tables IB1I and IB2I show the results from the fixed R e fits 
to our data, compared with the free radius fits and the best 
fits to the rest-frame optical images of K01. 



Object 


x 2 d 




jtte 


rnuc 
AB 


rhost 
AB 


PA 


a/b 


SGP2:36 


1.244 


1.6 


22.3 


21.5 


22.5 


136 


1.29 




1.246 


2.0 


22.9 


21.4 


22.5 


156 


1.28 




1.272 


5.1 


25.0 


21.3 


22.7 


152 


1.32 




1.300 


10.2 


26.5 


21.2 


22.7 


150 


1.39 


ff-band 




- 5 




21.3 


21.0 






SGP2:25 


1.295 


1.8 


23.2 


20.8 


23.1 


98 


1.21 




1.617 


2.0 


23.2 


20.9 


22.9 


120 


1.11 




1.635 


5.0 


25.7 


20.8 


23.4 


156 


1.24 




1.643 


10.0 


27.3 


20.7 


23.5 


166 


1.80 


ff-band 




< 10 




20.9 


21.2 






SGP2:11 


1.169 


1.0 


25.2 


20.6 


23.4 


168 


1.15 




1.264 


2.0 


23.4 


20.7 


23.0 


35 


1.14 




1.263 


5.0 


25.6 


20.6 


23.3 


34 


1.12 




1.264 


10.0 


27.1 


20.6 


23.3 


31 


1.11 


ff-band 




< 10 




20.3 


21.9 






SGP3:39 


1.324 


2.2 


22.8 


20.6 


22.5 


107 


1.13 




1.325 


2.0 


22.8 


20.6 


22.5 


141 


1.14 




1.344 


5.0 


25.2 


20.5 


22.9 


113 


1.25 




1.356 


10.0 


26.8 


20.4 


23.0 


108 


1.47 


ff-band 




< 10 




20.8 


21.1 






SGP4:39 


1.031 


1.5 


22.1 


20.7 


22.5 


81 


1.09 




1.031 


2.0 


22.9 


20.6 


22.6 


113 


1.11 




1.036 


5.1 


25.3 


20.5 


23.0 


115 


1.18 




1.041 


10.2 


27.0 


20.4 


23.1 


115 


1.25 


H-band 




< 10 




20.2 


22.9 






PKS1524 


1.198 


2.2 


22.0 


19.6 


21.5 


163 


1.31 




1.256 


2.0 


21.7 


19.6 


21.4 


163 


1.30 




1.266 


5.0 


24.1 


19.5 


21.8 


162 


1.46 




1.281 


10.0 


25.7 


19.4 


21.9 


162 


1.69 


iJ-band 




~ 5 




19.4 


20.7 






B2 2156 


1.335 


3.6 


22.3 


19.4 


20.7 


39 


1.10 




1.393 


2.0 


20.7 


19.4 


20.4 


28 


1.11 




1.347 


5.0 


23.1 


19.3 


20.8 


24 


1.12 




1.422 


10.0 


24.6 


19.1 


20.9 


20 


1.15 


ff-band 




12.1 




19.3 


19.2 




1.77 


PKS2204 


1.300 


1.6 


21.3 


20.8 


21.6 


137 


1.48 




1.364 


2.0 


22.0 


20.8 


21.7 


110 


1.42 




1.400 


5.0 


24.2 


20.5 


21.9 


138 


1.78 




1.431 


10.0 


25.8 


20.3 


22.0 


136 


2.31 


iJ-band 




- 5 




19.9 


22.0 






4C45.51 


1.767 


2.2 


22.0 


20.6 


21.5 


78 


1.13 




1.829 


2.0 


21.8 


20.7 


21.5 


76 


1.10 




1.859 


5.0 


24.0 


20.5 


21.7 


86 


1.10 




1.906 


10.0 


25.6 


20.4 


21.7 


91 


1.12 


H-band 




16.5 




18.8 


19.2 




1.23 
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